Scaling behavior of developing and decaying networks 



(N 
> 



O 

o 
o 



o 



S.N. Dorogovtseyi'^^* and J.F.F. Mendes^-t 
^ Departamento de Fisica and Centra de Fisica do Porto, Faculdade de Ciencias, Universidade do Porto 
Rua do Campo Alegre 687, 4169-007 Porto, Portugal 
^ A.F. loffe Physico-Technical Institute, 194021 St. Petersburg, Russia 

We find that a wide class of developing and decaying networks has scaling properties similar to 
those that were recently observed by Barabasi and Albert in the particular case of growing networks. 
The networks considered here evolve according to the following rules: (i) Each instant a new site 
is added, the probability of its connection to old sites is proportional to their connectivities, (ii) In 
addition, (a) new links between some old sites appear with probability proportional to the product 
; ^ - of their connectivities or (b) some links between old sites are removed with equal probability. 
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I Recently observed scaling behaviour of a number of networks [1 — 14] again sharply increased the interest of these 
exciting objects widely studied for a long time [15 — 20]. In fact the scaling properties were found only in a few of a 

■ great number of growing networks PJ4H7|,P^,?] but one of them - Web - is so significant for everybody that the topic 
^> . turned to be really hot. 

The simplest model of a scale- free growing network was proposed by Barabasi and Albert 0,^. In this model, each 
new site is connected with some old site with probability proportional to its connectivity fc, i.e. to the number of 
connections with this site. Then the distribution of the connectivities in the large network has a power-law dependence 
P{k) oc k~"' with the exponent 7 = 3 |@,||. In fact, such a network is self-organized into a scale-free structure. A full 
' form of the distribution of the connectivities and some other related properties of the model were calculated exactly 
iT) 10 . Introduction of aging of sites proportional to t^" , where t is the age of a site, does not change scaling properties 
crucially for a < 1, but scaling breaks at higher values of the aging exponent [ p3| . Several examples of real networks 
^ ■ with aging of sites are described in The simplicity of the Barabasi-Albert's (BA) model makes it a convenient 
object to study evolution of networks. 

Nevertheless, the BA model describes only a particular type of evolving networks. Of course, reality is much reacher. 
In real networks, e.g. in Internet, links are not only added but may break from time to time. That certainly changes 
' the structure of such networks. Note also, that new links between old sites may appear in a different way than links 
between new and old sites, according to different rules. Therefore, it is tempting to find out whether the observed 
I I behaviour is usual for a vast variety of networks or applies only to a very restricted number of invented objects. With 
''O ' that purpose, we extend a set of models of evolving networks starting from the BA model. 
Ch , In the present letter, we propose models of developing and decaying networks with undirected links which show 
scaling behaviour. We consider structures which evolve due to the following reasons. First, they grow like in the 
BA model, i.e. in each instant one new site is added and is connected with an old site by an undirected link with a 
probability proportional to its connectivity k (one may check that the general results ~ the existence of the scaling 
k>( , and the values of the scaling exponents - do not depend on the number of the connections with a new site). In 
j_j ■ addition, we introduce a new parallel component of the evolution - the permanent addition of new undirected links 

■ between old sites or, on the contrary, the permanent removal of some old links. We consider two different cases, 
(a) A developing network: Each instant, new c links are added between unconnected pairs of old sites i and j with 
probability proportional to the product of their connectivities kikj. (c may be also non integer. For that one can 
introduce probability of addition of a link.) c > 0. (b) A decaying structure (in fact, it is a set of clusters): Each 
instant, some links between old sites are removed with equal probability. In this case c < 0. 

Note that both processes - the addition of new sites with new links and the addition of new links between old sites 
(or removing old links) proceed in parallel, so the resulting structures differ from the original BA model all the time. 

We study the following one-site quantities of the structures: the total distribution of connectivities at long times, 
P(fc), and the average connectivity of a site of an age s at long time t, k{s, t), and their scaling exponents P{k) oc 
and k{s,t) oc (s/t)-'^. 

Below, we demonstrate both analytically and by simulation that the introduced evolving networks show scaling 
behaviour in a wide range of values of c. Nevertheless, while both P{k) and k{s,t), for the developing networks, are 
power-law functions for all c > 0, only k{s, t) demonstrates the power-law behaviour in the whole range — 1 < c < 
for the decaying structures. In this case, the power-law dependence of the distribution P{k) is observed only close to 
c = 0. 

In order to study scaling properties of the evolving networks we performed numerical simulations according to the 
introduced above rules. Each instant, we add one new site with one link and, in addition, may remove some of the 
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old links (decaying network) or, on the contrary, may add some new links between unconnected directly old sites 
(developing network) with the relative rate c. Therefore, to study even one-site properties of the structures, one has 
to keep in memory information about all connections among them. We performed simulations with a total number 
of sites (e. g. time) t = 1000 with 100000 averages for decaying networks, — 1 < c < 0, and t = 10000 with 10000 
averages for developing networks, c > 0. In Fig. 1, we present the dependences of the average connectivity k{s,t) on 
the number of a site s at different values of c for both structures, i.e. for the decaying network. Fig. 1 (a), and for 
the developing one. Fig. 1 (b). 

For both models, in the whole range of c, k{s,t) oc {s/t)~^. The change of the sign of the exponent (3 in the 
developing network at c = —1/2 was unexpected (compare with the behaviour of f3 vs. an aging exponent in networks 
with aging of sites (23|), see Fig. 1 (a) and the dependence /3(c) in Fig. 2. At this point, the average connectivity 
turns to be independent of the site age, k{s,t) = 1. We studied also the distribution P{k). It behaves as fc""*" for all 
c > for the developing network but, for the decaying network, the power-law dependence is found only in a narrow 
region of c near zero (see Fig. 2). The range of the values of k for which we observe such behaviour diminishes with 
decrease of c and then disappears. Note that the finite size effects are strong in this region (c < 0). 

We studied also the models in which, unlike the structures considered above, old links between sites are permanently 
removed with probability proportional to the product of the connectivities of the sites, or new links between old sites 
are permanently added with equal probability. The simulation demonstrates that the scaling breaks in both cases. 

Let us describe the obtained results analytically. We start from the case of the developing network. One may use 
the simple continuous approach [p|j23| that gives exact results for the scaling exponents as it was demonstrated in 
pl[ . Since one site is added per unit of time, then the total number of sites is t and each site is labeled by the time of 
its birth s < t. Then the equation for the average connectivity of the site s at time i, fc(s, i), in the continuous limit 
may be written in the following form: 
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(we wrote also the boundary condition - one link connected with a new site is added at each time step). The first 
term is the same as in the BA model, and the second one describes the increase of the connectivity due to the addition 
of new links between old sites with probability proportional to the product of connectivities of the connected sites. 
Note that the new links between the old sites can appear only if there is still no links between them. Hence, to write 
the last term in the present form, we make a strong assumption: we assume that the effect of multiplying of links is 
not essential at long times. We checked the validity of this non obvious assumption by simulation. 
Eq. (nl) is simplified at long times: 
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Applying ds to Eq. (||) one gets 
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dsk{s,t) - k{t,t) = 1 + 2c, 



so we obtain the obvious relation Jg dsk{s, t) ~ 2(1 + c)t. Now Eq. @) is of the following simple form: 

dk{s,t) _ l + 2c k{s,t) 
dt ~ 2(1 -he) t ■ 

It solution is k{s,t) = (s/t)^^ with the exponent 

l + 2c 
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To obtain the exponent of the distribution of connectivities, 7, one uses the general relation between the scaling 
exponents of growing networks |23,El| : 
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/?(7 - 1) = 1 
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that was obtained on the assumption that both fc(s, t) and P{k) show scahng behaviour. From our simulation, we 
know that this condition is fulfiUed for c > 0. Therefore, 
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l + 2c' 
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so we get both scahng exponents for the developing network. 

Let us consider now the decaying network. Again we apply the continuous approach. In fact, the removal of old 
links with equal probability seems to be equivalent to the decrease of the connectivities of old sites with probability 
proportional to their particular values, so Eq. (||) may also be applicable to this case. Nevertheless, one should 
account for the fact that only existing links may be removed. Therefore, we prefer to make the calculations more 
thoroughly. 

One introduces the average number of links between the sites s and s' at time t, n{s,s',t), where < s < s' < t. 
The average connectivity may be expressed in terms of this quantity: 



k{s,t) = / dun{u,s,t)+ / dwn{s,w,t) 
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The set of equations for n(s, s', t) is 
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(Note that c is negative now!) The first equality of Eq. (g) describes the links added to the network together with 
new sites as in the BA model. We again set the number of links connected with each new site to be unit. Applying 
ds to this equality we get 



/" 

^0 



ds n{s, s' , s') = 1 . 



(10) 



The second equality of Eq. (ph shows how n{s,s',t) changes due to the removing of links between the old sites. 



Application of ds ds' to this equality leads to the other obvious relation, ds ds' n{s, s' , t) = (1 + c)t 
Let us search the solution of Eq. (O) in the scaling form 
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and 



One sees that 
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and 



(l + c). (15) 



The solution of Eq. (n2h may be found in the form: 



where B,a,b are constants. Inserting Eq. ( p^ into Eq. ( p^ with account for Eqs. ( p^ and (^5|) we obtain the 
exponents a = b = 1 — 1/[2(1 + c)]. Eq. (|l^) gives B — I — a = 1/[2(1 + c)]. Substitution of Eq. (|16|) with account 
for Eq. ( Pt|) into Eq. (^) leads to the expression (3 = a = 1 — 1/[2{1 + c)] that is exactly the same as in the previous 
case, see Eq. (||). Note that now it is possible to use the relation between the scaling exponents, Eq. (^), only in the 
region near c = 0, since only in this region we observed the scaling behaviour of P{k). For these values of c we again 
get the old expression Eq. (j^) for the exponent 7. 

In Fig. 2, we plot the analytically obtained dependences (3 and 7 vs. c together with the results obtained from 
the simulation. One may see that the correspondence between the simulation and the theory is quantitative. When 
c changes from — 1 to 0, /? increases from —00 to 1/2 passing zero at c~ —1/2. Subsequent increase of c to 00 leads 
to growth of /? up to 1 while 7 decreases from 3 to 2. The particular case c cx3, 7 = 2, corresponds to the situation 
when the network evolves only due to the addition of new links by the defined above rules, that resembles the original 
small- world networks of Watts and Strogatz ^ . 

Our results show that the permanent removing of links leads to a more essential change of the structure of a network 
than the addition of them. What is the reason for that? One may see that the decaying structure under consideration 
is, in fact, a changing set of disconnected clusters. Because of finite size effects, we failed to find the position of the 
percolation threshold that may be defined for networks |J2|. Nevertheless, we see that, at high enough rates of link 
removal, large clusters are certainly absent, and the appearing structures indeed have to demonstrate quite different 
properties than the networks with c > 0. We failed also to find any peculiarity in the distribution of clusters in the 
point of the scaling break, c = —1/2. 

In summary, we have introduced a new parallel component of the evolution of growing networks. In addition to 
new links connecting new sites and old ones, links between old sites may appear or break with the relative rate c. We 
have demonstrated that addition of this component to a scale-free network does not break the scaling behaviour in a 
wide range of the rate c. 

The following questions remain open. What is the meaning of the point c = — 1/2 in which the exponent (3 changes 
sign? Is there any peculiar value of negative c below which the power-law behavior of P{k) breaks? 

SND thanks PRAXIS XXI (Portugal) for a research grant PRAXIS XXI/BCC/16418/98. JFFM was partially 
supported by the projects PRAXIS/2/2. l/FIS/299/94. We also thank M.C. Marques for reading the manuscript and 
A.V. Goltsev, Yu.G. Pogorelov and A.N. Samukhin for many useful discussions. 
* Electronic address: sdorogov@fc.up.pt 
^ Electronic address: jfmendes@fc.up.pt 



[1] J. Lahererre and D. Sornette, Eur. Phys. J. B 2, 525 (1998). 

[2] S. Redner, Eur. Phys. J. B 4, 131 (1998). 

[3] D. J. Watts and S. H. Strogatz, Nature 393, 440 (1998). 

[4] R. Albert, H. Jeong, and A.-L. Barabasi, Nature 401, 130 (1999). 

[5] B.A. Huberman and L.A. Adamic, Nature 401, 131 (1999). 



4 



B. A. Huberman, P. L. T. PiroUi, J. E. Pitkow, and R. J. Lukose, Science 280, 95 (1998). 

A.-L. Barabasi and R. Albert, Science 286, 509 (1999). 

A.-L. Barabasi, R. Albert and H. Jeong, Physica A 272, 173 (1999). 

D. J. Watts, Small Worlds (Princeton University Press, Princeton, New Jersey, 1999). 

J.B.M. Barthelemy and L.A.N. Amaral, Phys. Rev. Lett. 82, 3180 (1999); erratum: Phys. Rev. Lett. 82, 5180 (1999). 
M.E.J. Newman and D.J. Watts, Phys. Lett. A 263, 341 (1999). 



M.E.J. Newman and D.J. Watts, Phys. Rev. E 60, 7332 (1999); C. Moore and M.E.J. Newman, |cond-mat/0001393| , to be 
published. 

A. Barrat and M. Weigt, Eur. Phys. J. B 13 547 (2000). 

C. Tsallis and M.P. de Albuquerque, Eur. Phys. J. B 13, 777 (2000). 
P. Erdos and A. Renyi, Publ. Math. Inst. Hung. Acad. Sci. 5, 17 (1960). 

B. Bollobas, Random Graphs (Academic Press, London, 1985). 
G. Parisi, J. Phys. A 19, L675 (1986). 

B. Derrida and Y. Pomeau, Europhys. Lett., 1, 45 (1986). 
B. Derrida and D. Stauffer, Europhys. Lett. 2, 739 (1987). 
R. Monasson and R. Zecchina, Phys. Rev. Lett. 75, 2432 (1995). 

S.N. Dorogovtsev, J.F.F. Mendes, and A. N. Samukhin, con d-mat/0004..., to b e published. 



L.A.N. Amaral, A. Scala, M. Barthele my, and H.E. Stanl ey, |cond- mat/000145^ 



S.N. Dorogovtsev and J.F.F. Mendes, r:ond-mat/00014ie 



FIG. 1. Average connectivity vs. number s of the site for (a) decaying (— 1 < c < 0) and (b) evolving (c > 0) networks. 
One site per unit of time is added. Different curves correspond to different values of c. The network size is t = 10000 (a) and 
t = 1000 (b). 



FIG. 2. Exponents /3 of the average connectivity and 7 of the distribution of connectivities vs. c, i.e. vs. the rate of removal 
(c < 0) or addition (c > 0) links between old sites. Points are obtained from the simulations. The lines are found analytically 
[see Eqs. (^) and (j^)]. For the decaying network, the scaling behaviour of P{k) is observed only in narrow region of c close to 
zero. 
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